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Abstract 

In the model of local computation algorithms (LCAs), we aim to compute the queried part of 
the output by examining only a small (sublinear) portion of the input. Many recently developed 
LCAs on graph problems achieve time and space complexities with very low dependence on 
n, the number of vertices. Nonetheless, these complexities are generally at least exponential 
in d, the upper bound on the degree of the input graph. Instead, we consider the case where 
parameter d can be moderately dependent on n, and aim for complexities with subexponential 
dependence on d, while maintaining polylogarithmic dependence on n. We present: 

• a randomized LCA for computing maximal independent sets whose time and space com¬ 
plexities are quasi-polynomial in d and polylogarithmic in n; 

• for constant e > 0, a randomized LCA that provides a (1 — e)-approximation to maximum 
matching whose time and space complexities are polynomial in d and polylogarithmic in 

n. 
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1 Introduction 


In the face of massive data sets, classical algorithmic models, where the algorithm reads the entire 
input, performs a full computation, then reports the entire output, are rendered infeasible. To 
handle these data sets, the model of local computation algorithms (LCAs) has been proposed. 
As defined in [RTVXlT] , these algorithms compute the queried part of the output by examining 
only a small (sublinear) portion of the input. Let us consider the problem of finding a maximal 
independent set (MIS) as an example. The algorithm A is given access to the input graph G , then 
it is asked a question: “is vertex v in the MIS?” The algorithm then explores only a small portion 
of G , and answers “yes” or “no.” The set of vertices {v : A answers “yes” on u} must indeed 
form a valid MIS of G. LCAs have been constructed for many problems, including MIS, maximal 
matching, approximate maximum matching, vertex coloring, and hypergraph coloring ( [RTVX111 
IARVX121 IMRVX121 IMV131 IEMR14al IRV14j h In our paper, we study MIS and approximate 
maximum matching; these are fundamental graph problems, well-studied in many frameworks, and 
moreover, tools and results for these problems have proven to be useful as building blocks for more 
sophisticated and specialized problems in the field. 

The LCA framework is motivated by the circumstances where we focus on computing a small, 
specified portion of the output. This key characteristic of LCAs generalizes many other models from 
various contexts. For instance, LCAs may take the form of local filters and reconstructors [CS061 
IACCL081 IBra08l IKPS081 ISS101 IJR131 ICGR131 ILRR14j . Important applications include locally 
decodable codes (e.g., [STV99]), local decompression [MSZ05 , D LRR13 . and locally computable 
decisions for online algorithms and mechanism design |MRVX12l IHMV14] , There are a number 
of works on related models of local computation (e.g., [ACL061 lBBC + 12l IST131IOZ14] ) as well as 
lower bounds for the LCA framework and other models (e.g., |GS14| h 

Many recently developed LCAs on graph problems achieve time and space complexities with 
very low dependence on n, the number of vertices. Nonetheless, these complexities are at least 
exponential in d, the upper bound on the degree of the input graph. While these papers often 
consider d to be a constant, the large dependence on d may forbid practical uses of these algorithms. 
In this work we consider the case where the parameter d can be moderately dependent on n, and 
provide LCAs for complexities that have quasi-polynomial and even polynomial dependence on d , 
while maintaining polylogarithmic dependence on n. As noted in [MRVX12] . whether there exist 
LCAs with polynomial dependence on d for these problems is an interesting open question. Our 
paper answers this question for the approximate maximum matching problem in the affirmative, 
and aims at providing techniques useful towards resolving other problems. 

1.1 Related Work 

Many useful techniques for designing LCAs originate from algorithms for approximating the solution 
size in sublinear time. For example, if we wish to approximate the size of the minimum vertex cover 
(VC), we may sample a number of vertices and check whether each of them belongs to the minimum 
VC or not. The main difference, however, is that an LCA must be able to compute the answer to 
every query, while an approximation algorithm is not required to produce a consistent answer for 
every sample, and may use other properties of the problem to infer its answer. This paper makes 
use of a number of common techniques from these approximation algorithms. 

In our work we build on the Parnas-Ron reduction, proposed in their paper on approximating 
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the size of a minimum VC |prH 7| . This reduction turns a fc-round distributed algorithm into an 
LCA by examining all vertices at distance up to k from the queried vertex, then simulating the 
computation done by the distributed algorithm, invoking d°^ queries to the input graph in total. 
Using this technique, they obtain a 2-approximation for the size of a minimum VC (with en additive 
error) with query complexity d°^ og 0/ e )) ^ an( t a c-approximation for c > 2 using d°( lo s rf / e ) queries. 
Marko and Ron later improve this result to a (2 + (^-approximation with query complexity dP ^° grf ) 
|MR09| . Distributed algorithms for the MIS problem require more rounds, and consequently, a 
similar reduction only yields an LCA with query and time complexities d°^ dlogd ^ log n in jRTVXl 1| . 

Another powerful technique for bounding the query and time complexities is the query tree 
method from the Nguyen-Onak algorithm |NQ08| . This method aims to convert global algorithms 
that operate on the entire input into local algorithms that adaptively make queries when a new 
piece of information is needed. To illustrate this approach, let us consider the MIS problem as 
an example. Recall the sequential greedy algorithm where we maintain an initially empty set /, 
then iterate through the vertex set in some order, adding each vertex to I if it does not have a 
neighbor in I. Each vertex v will be in the resulting MIS if and only if none of v's neighbors 
that precede v in our order is already in the MIS. Prom this observation, we may create a local 
simulation of this algorithm as follows: to determine whether v is in the MIS, we make recursive 
queries to the preceding neighbors of v and check whether any of them is in the MIS. As we only 
query preceding neighbors, the structure of our recursive queries form a query tree. Nguyen and 
Onak apply this approach on a random order of vertices so that the size of the query tree, which 
determines the time and query complexities, can be probabilistically bounded. This method is used 
in [ARVX121 IMRVX121 IMV131 IHMV141 IRV14] , giving query complexities with polylogarithmic 
dependence on n for various problems. Unfortunately, the expected query tree size is exponential 
in d , which is considered constant in these papers. For certain problems, a slight modification of 
the Nguyen-Onak algorithm reduces the expected query tree size to 0(d) [YYI121 IORRR12] .This 
gives algorithms with poly(d) query complexity for approximating the sizes of maximum matching 
and minimum VC with multiplicative and additive errors. We build on these results in order to 
obtain LCAs whose query and time complexities are polynomial in d. 

Recently, a new method for bounding the query tree sizes using graph orientation is given in 
[EMR14a] based on graph coloring, which improves upon LCAs for several graph problems. They 
reduce the query complexity of their algorithm for the MIS problem to d°0 2 lo § rf ) log* n, giving the 
lowest dependence on n currently known, as well as a new direction for developing deterministic 
LCAs. This approach can also be extended back to improve distributed algorithms for certain cases 
[EMR.14bj . 

While all of these LCAs have complexities with exponential dependence on d for the problems 
studied in this paper, there has been no significant lower bound. To the best of our knowledge, 
the only lower bound is of Q(d), which can be derived from the lower bound for approximation 
algorithms for the minimum VC problem, given by Parnas and Ron |PR,07| . 

1.2 Our Contribution and Approaches 

This paper addresses the maximal independent set problem and the approximate maximum match¬ 
ing problem. The comparison between our results and other approaches is given in table [IJ Our 
paper provides the first LCAs whose complexities are both quasi-polynomial and polynomial in 
d and polylogarithmic in n for these problems, respectively. More concretely, when d is non- 
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Problem 

Citation 

Type 

Time 

Space 

MIS 

[RTVXll 

randomized 

20(dlog 2 d) \ 0 g n 

0(n) 

ARVX12 

randomized 

2 0(d\og 2 d) log 3 n 

2 0(dio g 2 d) log 2 n 

|EMR 14a| 

deterministic 

2 0(d 2 log 2 d) log * n f 

none 

[RVT4] 

randomized 

2°( d ) log 2 n 

2°( <i ) log n log log n 

2°( d ) log n log log n 

2°( d ) log 2 n 

this paper 

randomized 

2 *^ (log 3 d ) log 3 n 

2 0(log a d) log 2 n 

Approximate 

Maximum 

Matching 

MY 13 

randomized 

0( log 4 n)* 

0(log 3 n)* 

|EMR 14a,| 

deterministic 

2P°M rf ) poly (log* n )t 

none 

this paper 

randomized 

poly{d, log n} 

poly{d, log n} 


Tabic 1: The summary of complexities of various LCAs. For the approximate maximum matching problem, 
e is assumed to be constant, j indicates query complexity, when time complexity is not explicitly given in 
the paper. | indicates hidden dependence on d , which is at least 2 °C) but not explicitly known. 

constant, previously known LCAs have complexities with polylogarithmic dependence on n only 
when d = O(loglogn). Our LCAs maintain this dependence even when d = exp(0((log logn) 1//3 )) 
for the MIS problem and d = poly(logn) for the approximate maximum matching problem. 

1.2.1 Maximal Independent Set 

We provide an LCA for computing a MIS whose time and query complexities are quasi-polynomial 
in d. We construct a two-phase LCA similar to that of [RTVXlTj . which is based on Beck’s 
algorithmic approach to Lovasz local lemma |Bec91| . In the first phase, we find a large partial 
solution that breaks the original graph into small connected components. The LCA for the first 
phase is obtained by applying the Parnas-Ron reduction on distributed algorithms. The distributed 
algorithm for the first phase of the MIS problem in RTVX11J requires O(dlogd) rounds to make 
such guarantee. Ours are designed based on recent ideas from [BEPS12j so that 0(poly(log d)) 
rounds suffice. As a result, after applying the Parnas-Ron reduction, the query complexity on the 
first phase is still subexponential in d. Then, in the second phase, we explore each component and 
solve our problems deterministically; the complexities of this phase are bounded by the component 
sizes. By employing a technique from [ARVX12] . we reduce the amount of space required by our 
LCA so that it has roughly the same asymptotic bound as its time and query complexities. It is 
worth mentioning that our LCA for the MIS problem may be extended to handle other problems 
with reductions to MIS, such as maximal matching or (d + l)-coloring, while maintaining similar 
asymptotic complexities. 

1.2.2 Approximate Maximum Matching 

We provide an LCA for computing a (1 — e)-approximate maximum matching whose time and query 
complexities are polynomial in d. Our algorithm locally simulates the global algorithm based on 
Hopcroft and Karp’s lemma |HK73| . This global algorithm begins with an empty matching, then 
for 0(l/e) iterations, augments the maintained matching with a maximal set of vertex-disjoint 
augmenting paths of increasing lengths. Yoshida et al. show that this global algorithm has an 
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efficient local simulation in expectation on the queries and the random tapes |YYI12| . We derive 
from their analysis that, on most random tapes, this simulation induces small query trees on most 
vertices. From this observation, we construct an efficient two-phase LCA as follows. In the first 
phase, we repeatedly check random tapes until we find a good tape such that the query trees for 
most vertices are small. This test can be performed by approximating the number of vertices whose 
query trees are significantly larger than the expected size through random sampling. In the second 
phase, we simulate the aforementioned algorithm using the acquired random tape. On queries for 
which the query trees are significantly large, we stop the computation and report that those edges 
do not belong to the matching. Our good tape from the first phase limits the number of such edges, 
allowing us to acquire the desired approximation with high probability. 


2 Preliminaries 

2.1 Graphs 

The input graph G = ( V,E ) is a simple undirected graph with |Y| = n vertices and a bound on 
the degree d , which is allowed to be dependent on n. Both parameters n and d are known to the 
algorithm. Let d denote the average degree of the graph. Each vertex v G V is represented as a 
unique positive ID from [n] = {1,... ,n}. For v € V, let deg G (u) denote the degree of v, r G (u) 
denote the set of neighbors of v. and r G (u) = r G (F)U{u}. For U C V, define r G (I7) = U ue c/r G (u). 
The subscript G may be omitted when it is clear from the context. 

We assume that the input graph G is given through an adjacency list oracle O g which answers 
neighbor queries: given a vertex v £ V and an index i £ [d], the z th neighbor of v is returned if 
i < deg(u); otherwise, Y is returned. For simplicity, we will also allow a degree query which returns 
deg(u) when v is given; this can be simulated via a binary-search on O(logd) neighbor queries. 

An independent set I is a set of vertices such that no two vertices in / are adjacent. An 
independent set I is a maximal independent set if no other vertex can be added to I without 
violating this condition. 

A matching M is a set of edges such that no two distinct edges in M share a common endpoint. 
A matching is a maximal matching if no other edge can be added to M out violating this condition. 
Let V(M) denote the set of matched vertices, and \M\ denote the size of the matching, defined to 
be the number of edges in M. A maximum matching is a matching of maximum size. 

2.2 Local Computation Algorithms 

We adopt the definition of local computation algorithms from [RTVXllj . in the context of graph 
computation problems given an access to the adjacency list oracle O g . 

Definition 1 A local computation algorithm A for a computation problem is a (randomized) al¬ 
gorithm with the following properties. A is given access to the adjacency list oracle O g for the 
input graph G, a tape of random bits, and local read-write computation memory. When given an 
input (query) x, A must compute an answer for x. This answer must only depend on x, G, and 
the random bits. The answers given by A to all possible queries must be consistent; namely, all 
answers must constitute some valid solution to the computation problem. 
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The complexities of an LCA A can be measured in various different aspects, as follows. 

• The query complexity is the maximum number of queries that A makes to O g in order to 
compute an answer (to the computation problem) for any single query. 

• The time complexity is the maximum amount of time that A requires to compute an answer 
to any single query. We assume that each query to O g takes a constant amount of time. 

• The space complexity is total size of the random tape and the local computation memory 
used by A over all queries. 

• The success probability is the probability that A consistently answers all queries. 

In this paper, we refer to the time and query complexities rather exchangeably: while the time 
complexity may be much larger than the query complexity in certain cases, for all LCAs considered 
here, the time complexities are only roughly a factor of O(logn) larger than the query complexities. 
The space complexity of our LCAs are dominated by the size of the random tape, so we often refer to 
the space complexity as seed length instead. As for the success probability, we consider randomized 
LCAs that succeed with high probability, that is, the success probability can be amplified to reach 
1 — n~ c for any positive constant c without asymptotically increasing other complexities. 

2.3 Parnas-Ron Reduction 

Some of our algorithms apply the reduction from distributed algorithms to LCAs proposed by 
Parnas and Ron (PR07| . This reduction was originally created as a subroutine for approximation 
algorithms. Suppose that when a distributed algorithm A is executed on a graph G with degree 
bounded by d, each processor (vertex) v is able to compute some function / within k communication 
rounds under the LOCAL modelQ Then this function / must only depend on the subgraph of G 
induced by vertices at distance at most k from v. We can then create an LCA A' that computes / 
by simulating A. Namely, A' first queries the oracle to learn the structure of this subgraph, then 
makes the same decision on this subgraph as A would have done. In total, this reduction requires 
d°( fc ) queries to O g . 

2.4 Construction of Random Bits and Orderings 

While our LCAs rely on random bits and orderings, we do not require all bits or orderings to be 
truly random: LCAs tolerate some dependence or bias, as they only access a small portion of such 
random instance in each query. We now provide some definitions and theorems we will use to 
construct our LCAs. 

2.4.1 Random Bits 

We will use the following construction to generate k- wise independent random bits from the seed 
(truly random bits) given on the random tape. 

Theorem 1 ([ ABI86 ]) For 1 < k < m, there is a construction ofk-wise independent random bits 
x\,... ,x m with seed length O(fclogm). Furthermore, for 1 < i < m, each Xi can be computed in 
space 0{k\ogm). 

1 In the COCAC model, we optimize the number of communication rounds without limiting message size; in each 
round, each vertex may send an arbitrarily large message to each of its neighbors. 
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Note here that each random bit generated from this construction is either 0 or 1 with equal prob¬ 
ability. Nonetheless, for any positive integer q, we may generate a random bit that is 1 with 
probability exactly 1/q using O(logq) such truly random bits. 

2.4.2 Random Orderings 

For n > 1, let S n denote the set of all permutations on [n]. Some of our LCAs make use of random 
permutations of the vertex set. Generating such uniformly random permutations requires f l{n log to) 
truly random bits. Nonetheless, we apply the method from |AIWX12| to construct good random 
permutations for our LCAs. 

We generate our permutations by assigning a random value r(v) to each element v, and rank 
our elements according to these values. More formally, an ordering of [to] is an injective function 
r : [to] —> 1Z where 1Z is some totally ordered set. Let v±,... ,v n be the elements of [n] arranged 
according to their values mapped by r; that is, r(vi) < ••• < r(v n ). We call the permutation 
7T = (m,..., v n ) of [n] corresponding to this ordering the projection of r onto S n . We may refer to 
r(v) as the rank of v. In our construction, the size of the range of 1Z is polynomial in n. 

A random ordering V of [to] is a distribution over a family of orderings on [to]. For any integer 

2 < k < n, we say that a random ordering T> is k-wise independent if for any subset S C [n] of size 
k. the restriction of the projection onto S n of T> over S is uniform over all the k\ possible orderings 
among the k elements in S. A random ordering V is e-almost k-wise independent if there exists 
some k- wise independent random ordering T> such that the statistical distance between V and V 
is at most e. 

We shall use the following construction from Alon et al. [ARVX12] . 

Theorem 2 (i A RV X 12 4 Let n > 2 be an integer and let 2 < k < n. Then there is a construction 
of ( l/n 2 )-almost k-wise independent random ordering over [n] whose seed length is 0(fclog 2 n). 

3 Maximal Independent Set 

3.1 Overview 

Our algorithm consists of two phases. The first phase of our algorithm computes a large MIS, using 
a variation of Luby’s randomized distributed algorithm jLub86j . We begin with an initially empty 
independent set /, then repeatedly add more vertices to I. In each round, each vertex v tries to 
put itself into I with some probability. It succeeds if none of its neighbors also tries to do the same 
in that round; in this case, v is added to I, and is removed from the graph. 

By repeating this process with carefully chosen probabilities, we show that once the first phase 
terminates, the remaining graph contains no connected component of size larger than d 4 log n with 
high probability. This phase is converted into an LCA via the Parnas-Ron reduction. Lastly, in 
the second phase, we locally compute an MIS of the remaining graph by simply exploring the 
component containing the queried vertex. 
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3.2 Distributed Algorithm - Phase 1 


The goal of the first phase is to find an independent set / such that removing T + (I) breaks G into 
components of small sizes. We design our variation of Luby’s algorithm based on Beck’s algorithmic 
approach to Lovasz local lemma [BecQlj : this approach has been widely applied in many contexts 
(e.g., [Alo91l IKel93 , RTVX1H IBEPS12| L We design our algorithm based on the degree reduction 
idea from |BEPS12] Pl Our algorithm turns out to be very similar to the Weak-MIS algorithm 
from a recent paper by Chung et al. |CPS14j . We state their version, given in Algorithm [lj so that 
we may cite some of their results. 

Let us say that a vertex v is active if v (/ T + (J); otherwise v is inactive. As similarly observed 
in |Fel00] . applying a round of Luby’s algorithm with selection probability l/(d + 1) on a graph 
with maximum degree at most d makes each vertex of degree at least d/2 inactive with constant 
probability. To apply this observation, in each iteration, we first construct a graph G' of active 
vertices. Next, we apply Luby’s algorithm so that each vertex of degree at least d/2 becomes 
inactive with constant probability. We then remove the remaining high-degree vertices from G' 
(even if they may still be active). As the maximum degree of G is halved, we repeat this similar 
process for [log d] stages until G' becomes edgeless, where every vertex can be added to I. Since 
each vertex becomes high-degree with respect to the maximum degree of G' at some stage, each 
iteration gives every vertex a constant probability to become inactive. 


Algorithm 1 Chung et al.’s Weak-MIS algorithm 
1: procedure Weak-MIS(G, d) 

2: I i — 0 > begin with an empty independent set I 

3: for iteration i = 1,..., ci log d do > ci is a sufficiently large constant 

4: G' <— G[V \ T + (/)] > subgraph of G induced by active vertices 

5: for stage j = 1,..., [log d~\ do 

6: Vj •(— {w £ V(G') : deg G /(u) > d/ 2 J } > vertices with degree > half of current max. 

7: each v £ V(G') selects itself with probability pj = 1/(~=t + 1) > Luby’s algorithm 

8: if v is the only vertex in T G , (v) that selects itself then 

9: add v to I and remove T G ,(t>) from G' 

10: remove Vj from G' > remove high-degree vertices 

11: add V{G') to I > add lated vertices in G' to I 

12: return / 


Chung et al. use Weak-MIS to construct an independent set such that the probability that 
each vertex remains active is only 1/poly (d) [CPS14j . We cite the following useful lemma that 
captures the key idea explained earlier. The proof of this lemma is included in Appendix [Al 

Lemma 1 ( [CPS14] ) In Algorithm^ ifv £ Vj, thenv remains active after stage j with probability 
at most p for some constant p < 1. 

Observe that for v to remain active until the end of an iteration, it must be removed in step 10 
due to its high degree. (If v were removed in line 9 or line 11, then v would have become inactive 

2 Applying a similar reduction on the unmodified version gives an LCA with complexities 2°f log3 d+losdloslosn h 
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since either v or one of its neighbors is added to I.) Therefore, v must belong to one of the sets 
Vj. So, each vertex may remain active throughout the iteration with probability at most p. After 
ft (log d) iterations, the probability that each vertex remains active is only l/poly(d), as desired. 

Now we follow the analysis inspired by that of [BEPS12] to prove the guarantee on the maximum 
size of the remaining active components. Consider a set S CV such that distc(u, v) > 5 for every 
distinct u,v £ S. We say that S is active if every v £ S is active. As a generalization of the claim 
above, we show the following lemma. 

Lemma 2 Let S C V be such that dista(u,v) > 5 for every distinct u,v £ S, then S remains 
active until Weak-MIS terminates with probability at most c? S ^CI‘ S ’I). 

Proof: First let us consider an individual stage. Suppose that S is active at the beginning of 
this stage. By Lemma [TJ each vertex v £ S fi Vj remains active after this stage with probability 
at most p. Notice that for each round of Luby’s algorithm, whether v remains active or not only 
depends on the random choices of vertices within distance 2 from v. Since the vertices in S are at 
distance at least 5 away from one another, the events for all vertices in S are independent. Thus, 
the probability that S remains active after this stage is at most p\ SnV f. 

Now we consider an individual iteration. Suppose that S is active at the beginning of this 
iteration. By applying an inductive argument on each stage, the probability that S remains active 

Wl°g d] |cpjT/. I 

at the end of this iteration is at most p^i =1 1 Jl . Recall that for S to remain active after this 

iteration, every v £ S must belong to some set Vj. So, ^ ft Vj\ = I'S'l- Thus, S remains 

active with probability exp(—ft(|5|)). 

Lastly, we apply the inductive argument on each iteration to obtain the desired bound. H 
Now we are ready to apply Beck’s approach to prove the upper bound on the maximum size of 
the remaining active components ( |Bec91j . see also [RTVXlTj 1. 

Theorem 3 Weak-MIS(G, d) computes an independent set I of an input graph G within 0(log 2 d) 
communication rounds, such that the subgraph of G induced by active vertices contains no connected 
component of size larger than d 4 logn with probability at least 1 — l/poly(n). 

Proof: We provide a proof sketch of this approach. Let T be a tree embedded on the distance-5 
graph, defined as (V,{{u, v) : dist q(u,v) = 5}), such that distort, v) > 5 for every distinct u, v £ 
V(T). By Lemma El the probability that its vertex set S remains active is d~ n ( s \ where s = (S’l- 
It can be shown combinatorially that there are at most n(4d 5 ) s distinct trees of size s embedded 
on the distance-5 graph (for more details, see the proof of Lemma 4.6 in [RTVXlT] ). Therefore, 
the expected number of such trees whose vertex sets remain active is at most n(4d 5 ) s • d~^( s \ For 
s = logn, this quantity is bounded above by l/poly(n). By Markov’s inequality, all such sets S are 
inactive with probability at least 1 — 1/poly(n). 

Observe that any connected subgraph of G of size at least d 4 log n must contain some set S 
satisfying the aforementioned condition. Thus, the probability that any such large set remains 
active after Weak-MIS terminates is also bounded above by l/poly(n). ■ 

3.3 Constructing the LCA - Phase 1 

We now provide the Parnas-Ron reduction of Weak-MIS in Algorithm Q] into an LCA. Let us start 
with a single stage, given as procedure LC-MIS-Stage in Algorithm EJ The given parameters are 










the graph G' (via oracle access), the degree bound d, the queried vertex v, the iteration number i 
and the stage number j. Given a vertex v, this procedure returns one of the three states of v at 
the end of iteration i, stage j: 

• YES if v € I (so it is removed from G') 

• NO if v £ I and it is removed from G' 

• _L otherwise, indicating that v is still in G' 


Algorithm 2 LCA of Phase 1 for a single stage of Weak-MIS 
1: procedure LC-MIS-STAGE(0 G, ,d, u,i,j) 

2: if j = 0 then return _L o every vertex is initially in G' 

3: if LC-MIS-Stage^'- 17 , d, v, i, j — 1) ^ _L then 

4: return LC-MIS-Stage^^, d, v, i, j — 1) > v is already removed from G' 

5: for each u within distance 2 from v on G' do > check the status of every vertex near v 

6: if LC-MIS-STAGE(C? G, ,d,— 1) = _L then 

7: status(w) 4— removed 

8: else 

9: if B(u,i,j ) = 1 then status(w) 4— selected i> B(u,i,j ) = 1 means u chooses itself 

10: else status(u) 4— not selected 

11: if status(u) = selected AND \/u £ Tg"(u): status(u) / selected then 

12: return YES o v is added to I 

13: if status(u) = not selected then 

14: if 3u £ rc'(u): status(u) = selected AND \/w £ Tg/(u): status(?u) ^ selected then 

15: return NO > a neighbor u of v is added to I , so v is removed 

16: if |{u £ : status (it) / removed}] > d/2^ then 

17: return NO t> v is removed due to its high degree 

18: return _L > v remains in G' 


For simplicity, we assume that the local algorithm has access to random bits in the form of a 
publicly accessible function B : V x [ci logd] x [flog cT|] —> {0,1} such that B(v,i,j ) returns 1 with 
the selection probability pj (from Algorithm [T] line 7) and returns 0 otherwise. This function can 
be replaced by a constant number of memory accesses to the random tape, and we will explore 
how to reduce the required amount of random bits in Section 13.51 Note also that we are explicitly 
giving the oracle O g as a parameter rather than the actual graph G'. This is because we will 
eventually simulate oracles for other graphs, but we never concretely create such graphs during 
the entire computation. Observe that each call to LC-MIS-Stage on stage j invokes up to 0{d 2 ) 
calls on stage j — 1. Simulating all stages by invoking this function with j = [log d\ translates to 
^O(iogd) q uer j es the the base level j = 0. 

Next we give the LCA LC-MIS-Iteration for computing a single iteration in Algorithm [3j 
The parameters are similar to that of LC-MIS-Stage, but the return values are slightly different: 

• YES if v £ I (so it is inactive) 

• NO if v £ T+(I) (so it is inactive) 

• _L otherwise, indicating that v is still active 

In case v is still active by the end of iteration i — 1, we must simulate iteration i using LC-MIS- 
Stage. We must return YES if v £ I (line 7 of Algorithm [3]) ; which may occur in two cases. It 
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can be added to I in some stage, making LC-MIS-Stage returns YES (corresponding to line 9 of 
Algorithm |TJ). It may also remain in G' through all iterations, making LC-MIS-Stage return _L; 
v is then added to I because it is isolated in G' (line 11 of Algorithm [T]). 


Algorithm 3 LCA of Phase 1 for a single iteration of Weak-MIS 
1: procedure LC-MIS-Iteration(C g , d, v, i) 

2: if i = 0 then return _L > / is initially empty 

3: if LC-MIS-Iteration(C g , d, v, i - 1) / _L then 

4: return LC-MIS-Iteration(C g , d, v, i — 1) > v is already inactive 

5: O g ' •<— oracle for the subgraph induced by {u : LC-MIS-Iteration(C> g , d, u, i — 1) = _L} 

6: if LC-MIS-Stage((9 g \log d] ) / NO then 

7: return YES > v is added to / in some stage (YES) or at the end (A) 

8: else if 3u £ rc(v): LC-MIS-Stage(C g ,d,u,i, flog cT|) ^ NO then 

9: return NO > a neighbor of u of v is added to I 

10: return A > v is still active 


To implement this LCA, we must simulate an adjacency list oracle for G' using the given 
oracle for G, which can be done as follows. For a query on vertex v, we call LC-MIS- 
Iteration(0 , d, u, i — 1) on all vertices u at distance at most 2 away from v. This allows us 
to determine whether each neighbor of v is still active at the beginning of iteration i, as well as 
providing degree queries. We then modify the ordering of the remaining neighbors (by preserving 
the original ordering, for example) to consistently answer neighbor queries. That is, O g can be 
simulated through at most 0(d 2 ) function calls of the form LC-MIS-Iteration( 0 g , d, u, i — 1). 

Using this subroutine, the LCA LC-MIS-PhaseI for Weak-MIS can be written compactly as 
given in Algorithm 01 Via a similar inductive argument, we can show that each call to LC-MIS- 
PhaseI translates to d°^ og d ) = 2°( log ^ calls to the original oracle O g . The running time for the 
LCA is clearly given by the same bound. The memory usage is given by the amount of random bits 
of B, which is 0(n log 2 d). We summarize the behavior of this LCA through the following lemma. 

Lemma 3 LC-MIS-PhaseI is a local computation algorithm that computes the distributed 
Weak-MIS algorithm with time complexity 2°( log ^ and space complexity 0(n log 2 d). 


Algorithm 4 LCA of Phase 1 (Weak-MIS) 

1: procedure LC-MIS-Phase 1(0 g , d, v) 

2: return LC-MIS-Iteration((P g , d, v, c\ log d) 


3.4 Constructing the LCA - Phase 2 and the Full LCA 

Let G" be graph G induced by active vertices after Weak-MIS terminates. By Theorem [31 with 
high probability, G" contains no component of size exceeding d 4 logn. Therefore, to determine 
whether an active vertex v is in the MIS, we first apply a breadth-first search until all vertices 
in C(v), the component containing v, are reached. Then we compute an MIS deterministically 
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Algorithm 5 LCA of Phase 2 
1: procedure LC-Phase 2(O g ” , d, v) 

2: breadth-first search for d 4 log n steps on G" to find C v 

3: if \C V \ > d 4 logn then report ERROR i> only occurs with probability l/poly(n) 

4: deterministically compute an MIS Iq v of C v 

5: if v £ Ic v then return YES 

6: else return NO 


and consistently (by choosing the lexicographically first MIS, for instance). This procedure is 
summarized as LC-Phase2 in Algorithm [5j 

The algorithm only reports ERROR when the component size exceeds the bound from Theorem 
[31 Clearly, a call to LC-Phase2 makes at most poly(d)logn queries to O g in total. The lexi¬ 
cographically first maximal independent set of C v can be computed via a simple greedy algorithm 
with 0(d • ICy) time complexity. Overall, both time and query complexities are poly{d, log n}. 

Combining both phases, we obtain the LCA LC-MIS for computing an MIS as given in Algo¬ 
rithm [6l We now prove the following theorem. 


Algorithm 6 LCA for computing a maximal independent set 

1: 

procedure LC-MIS(0 g , d, v) 


2: 

if LC-MIS-PhaseI (O c ,d,v) / Y then 


3: 

return LC-MIS-Phase 1(0 g , d, v) 

c> v is already inactive 

4: 

O g " ■<— oracle for the subgraph induced by {u 

LC-MIS-Phase 1(0 g , d, u) = Y} 

5: 

return LC-MIS-PHASE2((I? G,, , d, v) 



Theorem 4 There exists a randomized local computation algorithm that computes a maximal in¬ 
dependent set of G with time complexity 2°( log log n and space complexity 0{n log 2 d). 

Proof: To obtain the time complexity, recall from Lemma [3] that each call to LC-PhaseI can be 
answered within 2 G( ' log ^ time using 2 G( ' log ^ queries to O g . Thus the adjacency list oracle O g ” 
can be simulated with the same complexities. The LCA for Phase 2 makes poly(d) log n queries 
to O g " , resulting in 2°( log log n total computation time and queries. The required amount of 
space is dominated by the random bits used in Phase 1. The algorithm only fails when it finds 
a component of size larger than d 4 logn, which may only occur with probability 1/poly(n) as 
guaranteed by Theorem [3[ I 

3.5 Reducing Space Usage 

In this section, we directly apply the approach from |ARVX12] to reduce the amount of random 
bits used by our LCA, thus proving the following theorem. 

Theorem 5 There exists a randomized local computation algorithm that computes a maximal in¬ 
dependent set of G with ra seed of length 2°( log ^ log 2 n and time complexity 2°( log (!> log 3 n. 
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Proof: Observe that the MIS algorithm in Theorem [3] constructed throughout this section does 
not require fully independent random bits in function B. Our algorithm can answer a query for any 
vertex v by exploring up to q = 2°( log ^ log n vertices in total. So, random bits used by vertices 
not explored by a query do not affect our answer. One bit is used by each vertex in each round 
of Phase 1, and thus we only create b = 0( log 2 d) random bits for each vertex. Therefore, out of 
0(n log 2 d) bits from function B, only q ■ b bits are relevant for each query. 

To generate random bits for function B, we will apply Theorem 12.4.21 As we may require bits 
of function B to be 1 with probability as low as 1 /{d + 1), we will need up to [logd] bits from 
this construction to obtain one bit for function B. So in our case, we have k = [logd] • q ■ b = 
2°(i° g d ) log n and m = [log d\ ■0(n log 2 d) = 0(n log 3 d). Thus the amount of space can be reduced 
to 0(k\ogm ) = 2°( log logn • log(nlog 3 d) = 2°( log log 2 n. Similarly, the required amount of 
time for computing each random bit becomes 2°^ og d ) log 2 n. H 

3.6 Other Remarks 

The proposed LCA for the maximal independent set problem can be used as a building block for 
other problems. For example, a (d+ l)-coloring can be computed by finding a maximal independent 
set of G x Kd +1 where K^+i is a clique of size d+ 1, resulting in an LCA within the same asymptotic 
complexities (see also |Lin92j ). 

We may apply our MIS algorithm to the line graph L(G) in order to compute a maximal 
matching on G within the same asymptotic complexities as Theorem [5j Nonetheless, Barenboim 
et al. provide a distributed routine for computing a large matching that yields a similar bound as 
Theorem [3] that requires only O(log c() rounds [BEPS12] . We apply the Parnas-Ron reduction on 
their algorithm as similarly done in Section[3]to obtain the following theorem, where the dependence 
on d is reduced from 2°9°g d ) to 2°9°s Additional details for this result can be found in 
Appendix [Bj 

Theorem 6 There exists a randomized local computation algorithm that computes a maximal 
matching of G with seed of length 2 0( - log d ^> log 2 n and time complexity 2°9°g d ) log 3 n. 

4 Approximate Maximum Matching 

4.1 Overview 

In this section, we aim to construct an LCA that provides a (1 — e)-approximation to maximum 
matching. To do so, we first address the simpler problem of computing a (1 — e)-approximation to 
MIS. We present this LCA because it allows us to explain our LCA for the approximate maximum 
matching problem in a much clearer manner, as they share both their main principles and the two- 
phase structure. Moreover, this LCA is useful as a subroutine for other problems. For example, an 
approximation to MIS of the line graph L(G) readily yields a (1/2 — e)-approximation to maximum 
matching of G. 
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4.2 LCA for Computing an Approximate Maximal Independent Set 

We now describe an LCA that provides a (1 — e)-approximation to MIS. More specifically, with 
high probability, the set of vertices I that our LCA answers YES must be a subset of some maximal 
independent set I satisfying |/| > (1 — e)\I\. Our algorithm is based on a local simulation of the 
classical greedy algorithm for computing an MIS. This greedy algorithm iterates over the vertex 
set according to some arbitrary order and adds a vertex to the constructed independent set if and 
only if none of its neighbors has been added previously. 

Algorithm [7] summarizes the local simulation of the greedy algorithm as suggested by Nguyen 
and Onak and further analyzed by Yoshida et al. [NO08 . [YYI12j . In addition to the adjacency 
list oracle O g and the queried vertex v, LS-MIS also receives as an input a permutation ir on the 
vertices. The set of vertices v such that LS-MIS(0 , vr, v) returns YES is the lexicographically first 
MIS according to tt, which is also the MIS outputted by the greedy algorithm when the vertices 
are iterated in this order. 

Algorithm 7 Local simulation of the greedy MIS algorithm 

1: procedure LS-MIS(0 G , tt, v) 

2: query 0 G for all neighbors of v 

3: let v \,..., Vk be the neighbors of v sorted according to tt 

4: for i = 1,..., k do 

5: if Vi precedes v in vr then 

6: compute LS-MIS(O g , 7T, Uj) 

7: if LS-MIS(0 g , 7 t, Vi) = YES then return NO 

8: return YES 


We now give an overview of the construction of our LCA. From the analysis of Yoshida et al. we 
obtain that in expectation over the queries and the random permutations, the query complexity 
of Algorithm [3 is polynomial in d, the degree bound [YYI12] . More formally, let R^(v) denote 
the number of (recursive) calls to LS-MIS during the evaluation of LS-MIS(0 g , tt, v). This result 
from Yoshida et al. can be succinctly stated as follows. 

Theorem 7 ([ YYI12 ]) For any graph G = (V, E ) with n vertices and rn edges, 

E?r eS n ,vev[Rn( v )\ — 1 + — • 

However, the query complexity of an LCA is determined by the maximum number of queries 
invoked in order to answer any single query, so this expected bound does not readily imply an 
efficient LCA. We may alter Algorithm [7] so that it computes an independent set that is relatively 
large but not necessarily maximal as follows. If a query on a vertex v recursively invokes too many 
other queries, then the simulation is terminated and the algorithm immediately returns NO, i.e., 
that v does not belong to our independent set. We show that Theorem 0 implies that for at least a 
constant fraction of the orderings, this modification yields a good approximation ratio. Nonetheless, 
an LCA must succeed with high probability rather than only a constant probability. We resolve 
this issue by repeatedly sampling permutations until we obtain a sufficiently good one. 

While generating a truly random permutation tt requires f l(n) random bits, we show that we 
may generate sufficiently good permutations by projecting random orderings from a distribution of 
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small entropy onto S n . Using the construction from Alon et al. as stated in Theorem [2 it follows 
that a small seed whose size is polylogarithmic in n suffices for our algorithm. From these outlined 
ideas, we are now ready to prove the following result. 

Theorem 8 There exists a randomized (1 — e) -approximation local computation algorithm for max¬ 
imal independent set with random seed of length 0((d 2 /e 2 ) log 2 n log log n) and query complexity 
0((d 4 /e 2 ) log 2 n log log n). 

Proof: Our algorithm builds on Algorithm [7] and consists of two phases as follows. The input of 
the algorithm is a random seed (tape) s and a queried vertex v £ V. In the first phase, using s 
and graph queries, we find a good permutation tt = ttg(s) over the vertex set. Note that this good 
permutation -k is dependent only on s and G, and therefore must be the identical for any query. 
Then in the second phase, we decide whether v belongs to our independent set using the altered 
local simulation that limits the number of invoked queries. 

Formally, we say that n £ S n is good if Pr v€ y [R//(v) > (\ < 7 where 7 = e/d, £ = 6 t/e and 
t = m/n + 1. Let 5 denote the error probability. We claim that, given access to D£, an £-wise 
independent ordering constructed from s, Algorithm [8] computes a good permutation tt with the 
desired success probability. This algorithm checks each constructed permutation 77 whether it is 
good by approximating the fraction p 7Ti of vertices whose induced query trees have size exceeding £. 
The set S of sampled vertices is sufficiently large that we obtain a good approximation p ni of p ni , 
and note that checking whether R~. (v) > £ can be accomplished by simply simulating Algorithm [7] 
until £ recursive calls are invoked. 


Algorithm 8 Phase 1 of the LCA for the approximate MIS algorithm (finding a good ordering) 
1 : procedure LC-AMIS-Phase 1(0 g , s) 

2: let S be a multi-set of 0(log(l/<5) log log(l/5)/7 2 ) vertices chosen uniformly at random 

3: for i = 1,..., 0(log(l/5)) do 

4: let 77 be the projection of an ordering independently drawn from D 1 

5: let p ni denote the fraction of vertices v € S for which Rff.(v) > £ 

6: if p 7Ti < 37/4 then 

7: return 7r = 77 

8 : report ERROR 


Now we show that with probability at least 1 — 5, Algorithm [8] finds a good tt. Let H 7T;V be 
the indicator variable such that H nv = 1 when R^(v) > £, and H 7rv = 0 otherwise. By Markov’s 
inequality and Theorem [7] it holds that 

Pryr e s n ,vev[Rn (v) > £] < 7/6 • 


That is, 


Evev[^TreS n [H-jr,v]\ < 7/6 • 


( 1 ) 


Next we aim to justify that we may obtain our random permutations by projecting orderings 
from Djl instead of drawing directly from S n . More specifically, we claim that for a fixed v, is 
identically distributed regardless of whether tt is drawn uniformly from S n or obtained by projecting 
an ordering r drawn from D^. To establish this claim, consider the evaluation of H W:V by a decision 
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tree Tc. v , where each inner vertex of 7 g,v represents a query to the rank r(u) identified with the 
queried vertex u € V. In order to evaluate H n:V we proceed down the tree Ta, v according to ir 
until we reach a leaf. After at most £ queries to 7 r, the value of the indicator H n ^ v is determined. 
Therefore, the depth of Tg,v is a t most l , and each leaf is associated with either 0 or 1. In other 
words, each leaf is identified with a sequence of at most t vertices and their corresponding ranks. 

Let si and s 2 be a pair of leaves such that: (1) s\ and S 2 are identified with the same sequence 
of vertices, v\,... ,Vk', and (2) the ranks of iq,..., Vk in si and s 2 induce the same permutation on 
v\,... ,Vk- Clearly the values corresponding to leaves si and s 2 are identical. Therefore, our claim 
follows from the definition of Awise independent ordering. Thus, from equation (JTJ) , we obtain 


Kev\EreD t [H*,v\] < 7/6 • ( 2 ) 

Next, we say that n 6 S n is very good if Pr„ e y \B^{v) > £\ < 7 / 2 . By Markov’s inequality and 
Equation ([2]) . an ordering r drawn from D 1 induces a very good permutation with probability at 
least 2/3. Therefore an ordering r drawn according to a distribution that is (l/n 2 )-almost Cwise 
independent random ordering over [n\ induces a very good permutation with probability at least 
2/3 — 1/n 2 . So by generating 0(log(l/<5)) permutations 77 ’s from such random ordering, with 
probability greater than 1 — 6/2, at least one of them must be very good. 

Let Ei denote the event that none of the selected 77 ’s is very good, so E\ occurs with probability 
smaller than 6/2. Let E 2 denote the event that there exists a selected 77 such that \p 7Ti —p ni | > r jn/ 4. 
By Chernoff’s bound and the union bound, E 2 occurs with probability at most 6/2. Given that E 2 
does not occur and that p ni < 3y/4, it follows that 77 is good; i.e., p n . < 7 . Therefore, given that 
both Ei and E 2 do not occur, Algorithm [ 8 ] indeed finds a good ordering. By the union bound, we 
obtain a good ordering 7 r with probability at least 1 — 6, as desired. 

Theorem [2] implies that a random seed of length 0((d/e) log 2 n) suffices in order to obtain a 
(l/n 2 )-almost Cwise independent random ordering over [n]. Since Algorithm [ 8 ] draws 0(log(l/<5)) 
such random orderings, overall a seed of length 0((d/e) log(l/<5) log 2 n) suffices. Additionally, we 
need 0((d 2 /e 2 ) log(l/<5) log log(l/<5) log n) more random bits to determine the set S of sampled 
vertices for evaluating p 7T . For Algorithm [ 8 ] to succeed with high probability, it suffices to substitute 
6 = l/poly(n), which yields the seed length claimed in the theorem statement. 

Finally, we now turn to describe the second phase of the algorithm. Let I n denote the maximal 
independent set greedily created by 7 r G S n , as defined by Algorithm [71 Consider the following 
LCA for a good permutation it. On query v, simulates the execution of LS-MIS(0 ,77 v) 
until it performs up to £ recursive calls to LS-MIS. If LS-MIS(C ,G , 77 v ) terminates within £ steps, 
£ 7 ,- returns the answer given by LS-MIS(C2 G , 77 v). Otherwise it simply returns NO. 

Let I' n = In fl V/ where V/ = {v : Bf//{v) < £}, and notice that C n answers YES precisely on 
V. Since ir is good, we have that V/ > (1 — q)n and so \I' n \ > \I n \ — qn. Since for any maximal 
independent set I it holds that |/| > n/d, we obtain \I' n \ > (1 — e)\I n \. In other words, C n computes 
a (1 — e)-approximate MIS, as required. I 

4.3 LCA for Computing an Approximate Maximum Matching 

Now we turn to describe an LCA that provides a (1 — e)-approximation to maximum matching. 
Our LCA locally simulates the global algorithm based on the following theorem by Hopcroft and 
Karp: 
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Lemma 4 \HK73\ Let M and M* be a matching and a maximum matching in G. If the shortest 
augmenting path with respect to M has length 2i — 1, then \M\ < (1 — l/i)\M*\. 

This global algorithm is used in the many contexts, including distributed computing and ap¬ 
proximation (e.g., [LPSP081INQ081 IYYH'2| 1. This algorithm can be summarized as follows. Let 
Mo be an empty matching, and let k = [1/e]. We repeat the following process for each i = 1 ,k. 
Let Pi denote the set of augmenting paths of length 2i — 1 with respect to Mj_i ■ We compute a 
maximal set Ai of vertex-disjoint augmenting paths, then augment these paths to M,._i to obtain 
Mi] that is, we set Mj = Mj_iAAj where A denotes the symmetric difference between two sets. 
As a result, Mk will be a (1 — e)-approximate maximum matching. 

The local simulation of this algorithm uses the formula M, = M*_i AAj to determine whether the 
queried edge e is in M in a recursive fashion. To determine Ai, consider the following observation. 
Let Hi = ( Pi,Ei ) be the graph whose vertex set corresponds to the augmenting paths, and (u,v) £ 
Ei if and only if the augmenting paths u and v share some vertex. Ai is a maximal set of vertex- 
disjoint augmenting paths; that is, Ai is an MIS of Hi. In order to compute locally whether a path 
is in Ai or not, Yoshida et al. apply Algorithm [7] on H,. They show that in expectation over the 
queries and the random permutations of paths in P\..... P^, the query complexity is polynomial 
in d. More formally, let n = (- 7 T 1 ,..., 7 r fc ) where ir l is a permutation of vertices (augmenting paths) 
in Pi, and let k (v) denote the query complexity of this local simulation. Yoshida et al. prove 
the following theorem. 

Theorem 9 ([ YYI12 ]) For any graph G = (V, E ) with n vertices and maximum degree d, 

er£,*M] < d“ 3 fc 0( ‘>, 

where tt = (-7T 1 ,..., 7r k ) and the expectation is taken over the uniform distribution over S , ip 1 i x ... x 
S \P k \ x 1^1- 

Based on this theorem, we convert the aforementioned local simulation into an LCA as similarly 
done in the previous section for the approximate MIS problem. We now prove the following theorem. 

Theorem 10 There is a randomized (1 — e)-approximation local computation algorithm for max¬ 
imum matching with random seed of length 0(d ek2+1 k 0 ^ log 3 n log log n) and query complexity 
0 (^ 6 fc 2 + 2 £.O(fc) i 0 g 2 n l 0 gl 0 g n ) where k = 0 (l/e). 

Proof: We build on the local algorithm of Yoshida et al. ( jYYI12j ). but here we double the 
number of iterations to k = [ 2 /e] so that Mk is a (1 — e/ 2 )-approximate maximum matching. 
The proof for this theorem follows the structure of the proof of Theorem [ 8 ] Our LCA contains 
two phases: the first phase finds a sequence of good permutations, then the second phase locally 
simulates the approximate maximum matching algorithm using those permutations. We say that 
a sequence of permutations It = {tt\, ... ,i r^.) £ S'|v 1 | x ... x S]y fc | is good if Pr ee p [Q?(e) > £\ < 7 
where 7 = e/( 2 d), I = 6 t/q and t = d 6fc2 k°^ k \ 

In the first phase, rather than individually generating a permutation 7 ^’ of It on each Pi, we create 
a single permutation 7 r over [m n} k\ where m Ht k = Yli=\ (t) — Yli =1 I P%\- I n terms of constructing 
random orderings, this simply implies extending the domain of our orderings to cover all augmenting 
paths from all k iterations. Then for each i £ [k\, the permutation ir k of paths Pi can be obtained 
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by restricting if to Pi (i.e., considering the relative order among paths in Pi). Since the algorithm 
queries if on at most i locations, we conclude that Theorem [9] holds for any if that is Ywise 
independent. Similarly to the proof of Theorem [TUI with probability at least 1 — l/poly(n), we 
shall find a good if by testing 0(logn) random orderings obtained via the construction of Alon 
et al. in Theorem [2j Therefore, to this end we need at most 0(Tlog 3 n + (d 2 /e 2 ) log 2 n log log n) 
random bits, as required. 

Now that we obtain a good sequence of permutations if from the first phase, we again perform 
the altered version of local algorithm that returns NO as soon as the simulation invokes too many 
queries. Let M.% denote the matching obtained by the algorithm of Yoshida et al. when executed 
with the permutations from if. Let ML = M^ n EL where EL = {e : Q~(v) < £}. Similarly to the 
proof of Theorem^ given if, we obtain an LCA that answers according to ML with query complexity 
i. If if is good, we have that EL > (1— r y)\E\ and so \ML\ > \M^\—j\E\ > (1 — e/2)|M*|— j\E\ where 
M* denotes a maximum matching. Since for any maximal matching M it holds that \M\ > \E\/{2d) 
we obtain \ML\ > (1 — e)|M*|, as required. ■ 
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A Proof of Lemma [T] 

Lemma 1 ([ CPS14 ]) In Algorithm^ ifv € Vj, thenv remains active after stage j with probability 
at most p for some constant p < 1. 

Proof: For each u G Tq>(v), let E u denote the event where u is the only vertex in L^, ({u, u}) that 
selects itself. Since the maximum degree in G' is at most d/ 2 J " —1 , then 

Pr[E u ] = Pj (l> Pj (l-pjf ^- 1 > g. 

Notice that E u is disjoint for each u G r G /(u). Since v G Vj, then deg G /(u) > d/2T Thus v 
becomes inactive with probability at least 

E Pr t^] ^ de SG' i v ) • (J) > ^2- 

«er G , (v) 

That is, the theorem holds with parameter p = 1 — l/4e 2 . ■ 

B Details of Theorem [6] 

Theorem 6 There exists a randomized local computation algorithm that computes a maximal 
matching of G with seed of length 2°^ og log 2 n and time complexity 2° dog log 3 n. 

In jII86j . Israeli and Itai proposed a randomized distributed algorithm which takes O(logn) 
rounds to compute a maximal matching. Similarly to the MIS problem, Barenboim et al. also 
create a variant of this algorithm that, within O(logd) rounds, finds a large matching that breaks 
the remaining graph into small connected components [BEPSl2j . Specifically, by running Algorithm 
1 the remaining graph satisfies the following lemma. 
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Algorithm 9 Barenboim et al.’s variant of Israeli and Itai’s algorithm for computing a partial 
matching (simplified) 

1: procedure Distributed-MM-Phase1(G, d) 

2: initialize matching M = 0 

3: for i = 1,... ,C 2 log d do > C 2 is a sufficiently large constant 

4: initialize directed graphs F\ = ( V , 0) and F 2 = ( V , 0) 

5: each vertex s chooses a neighbor t (if any) uniformly at random, 

then add ( s,t ) to E{F\) 

6: each vertex t with positive in-degree in F\ chooses a vertex s 6 {s' : ( s',t ) £ E(Pi)} 

with highest ID, then add (s,f) to E(F 2 ) 

7: each node v with positive degree in F 2 chooses a bit b(v ) as follows: 

if v has an outgoing edge but no incoming edge, b(v) = 0 
if v has an ingoing edge but no outcoming edge, b(v) = 1 
otherwise, chooses b(v ) £ {0,1} uniformly at random 
8: add every edge (s,t) £ E(F 2 ) such that b(s) = 0 and b(t) = 1 to M, 

and remove matched vertices from G 


Lemma 5 ( |BEPS12] I Distributed-MM-Phase1(G, d) computes a partial matching M of an 
input graph G within O(logd) communication rounds, such that the remaining graph contains no 
connected component of size larger than 0(d A \ogn) with probability at least 1 — l/poly(n). 

The proof of this lemma also makes use of Beck’s analysis, but contains a more complicated 
argument which shows that with constant probability, each vertex loses some constant fraction of 
its neighbors in every round. Thus, applying such matching subroutine for O(logd) rounds suffices 
to remove or isolate each vertex with probability 1 — 1/poly (d). We convert this lemma into a 
two-phase LCA in a similar fashion to obtain the desired LCA. 
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